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I review recent results regarding the nucleon charge and magnetization densities as well as the 
shape of the nucleon. First, some phenomenolgical considerations that show that the shape of 
the proton is not round are discussed. Then model independent results regarding the neutron and 
proton charge density, and the proton magnetization density are presented. Finally, I show how the 
spin-dependent densities that reveal the shape of the proton can be measured via their relation with 
transverse momentum distributions. The present work is made possible by the recent tremendous 
experimental progress made at Bates, Mainz and Jefferson Laboratory. 
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1. Phenomenology 

The electromagnetic form factors are matrix elements of the current operator, J ti {x), between 
nucleon states of different momentum: 

(p\XV^0)\pA)=u(p\X')^F l (Q 2 )+i—q a F 2 (Q 2 )y(p,X), (1.1) 

where q a = p' a — p a , Q 2 = —q 2 > 0, M is the nucleon mass and the light-cone helicities are A, A'. 
The Sachs form factors are: G E (Q 2 ) = F^Q 2 ) - £,F 2 (Q 2 ),G M (Q 2 ) = F y (Q 2 ) +F 2 (Q 2 ). The 
early expectation, based on a simple application of helicity conservation at very high values of the 
momentum transfer, was that the QF 2 (Q 2 )/2MF l ~ 2*=* _> ggi = const, which is also obtained 
from non-relativistic considerations. Thus the expectation was that, at sufficiently large momentum 
transfer, Ge/Gm would be flat and the QF 2 /F\ would fall. However, the reverse was true [jl|]. 

It is necessary to comment on the meaning of these form factors. In the Breit frame, with 
p = — p', Ge is the nucleon helicity flip matrix element of J°. However, any probability or density 
interpretation of Ge is spoiled by a non-zero value of Q 2 , no matter how small. The initial and final 
states have different momentum, and therefore relativistically have different wave functions. Any 
attempt to analytically correct for the total momentum by incorporating relativistic collections in a 
p /m 2 type of expansion is doomed by the presence of the very light current quark mass, m q . That 
is, at small values of Q , one finds 




(1.2) 



where the first term represents the expected effect depending on the square of the wave function 
and the term C represents the unknown boost correction. 

So to analyze form factors using a model, the model must be relativistic. We chose [§] to 
use a relativistic model using light front coordinates. These useful coordinates involve the use 
of a "time" x + = (ct + z)/V2 = (x° + x i )/\fl. The corresponding evolution operator is the not 
the Hamiltonian, p , but instead p~ = (P° — p 3 )/v / 2- The orthogonal spatial coordinate is x~ = 
(x Q —x 3 )/ y2. If one quantizes at x + = 0, then x~ = \[2z, and this why x~ is thought of as the 
spatial variable. The canonically conjugate momentum is given by p + = (p° + p 3 )/\/2- We note 
that p^ = p~x + + p + x~ — p b. The transverse coordinates perpendicular to the and 3 directions 
are denoted as b and p. Using these variables allows the separation of center of mass and relative 
position variables in a manner similar to that of the usual non-relativistic treatment. The key feature 
is that transverse boosts act like the non-relatistic boosts. 

We used these variables to formulate a Poincare invariant, relativistic constituent quark model 
in 1995 and predict the qualitative behavior of form factors measured five years later. Please see 
Figs. 10, 11 of Ref. j2|]. More detailed analysis ||3|]showed that the flat nature of ^ results from 
orbital angular momentum of the quarks inherent in the lower component of quark Dirac spinors. 

So we had a reasonable model, which includes quark orbital angular momentum. I was faced 
with the challenge of relating orbital angular momentum to a potential non-spherical shape of the 
proton. The notion that the proton might not be a sphere has its impetus in the discovery that the 
spins of quarks and anti-quarks account for only about 30% of the total angular momentum. It 
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seems natural to associate non-zero orbital angular momentum with a non-spherical shape, but the 
Wigner-Eckart theorem says that the proton has no quadrupole moment. In response, I introduced 
spin-dependent densities SDD, and later learned these are common in condensed matter physics. 
Spin-dependent density operators We interpret orbital angular momentum in terms of the shapes 
of the proton exhibited through the rest-frame ground-state matrix elements of spin-dependent 
density operators ph. The density operator is p (r) = 5(r — r,), where r, is the position operator 
of the i'th particle. For particles of spin 1/2 one can measure the combined probability that a particle 
is at a given position r and has a spin in an arbitrary, fixed direction specified by a unit vector n. 
The coordinate-space spin-dependent density SDD operator is p (r, n) = £ ; 5 (r — r,-) \ ( 1 + C; • n) . 

To understand the connection between the spin-dependent density and orbital angular mo- 
mentum, consider a first example of a single charged particle moving in a fixed, rotationally- 
invariant potential in an energy eigenstate \*¥i t i/2j} °f quantum numbers: I = l,j = 1/2, po- 
larized in the direction is and radial wave function R(r p ). The wave function can be written as 
( T p\ x Pi,lfo)=R(r p )<T-i p \s). The ordinary density. p(r) = (^i,!^ 1 5 (r - r p ) ^1,1/2^ > = R 2 {r), a 
spherically symmetric result because the effects of the Pauli spin operator square to unity. But the 
matrix element of the SDD is more interesting: 



The magnetic quantum defines an axis, s and the direction of vectors can be represented in terms 
of this axis: s • f = cos 6. Suppose n is either parallel or anti-parallel to the direction of the proton 
angular momentum vector s. Then p(r,n = s) =R 2 (r) cos 2 6, p(r,n = — s) = R 2 (r) sin 2 6, and the 
non-spherical shape is exhibited. The average of these two cases is a spherical shape. 

We computed momentum-space SDDs [Q] using our model|Q] model and obtained a variety of 
unusual shapes. See Figs. 2,3 of The logic is that using a model wave function, in rough agree- 
ment with data for form factors, leads to a non-spherical shape of the SDD. However, there was 
no direct connection between experiment and the shapes. We shall return to issue after discussing 
other model independent results for charge and magnetization densities. 

2. Model independent neutron charge density 

The neutron has no net charge, but the charge density need not vanish. So we can ask, "Is the 
central charge density negative or positive?". Long-standing existing answers are based on models 
For example,, the neutron can make a spontaneous quantum transition to a pn~ state []5]]. The light 
pion spreads out over a larger region of space than the proton to cause a negative charge density 
at the edge of the neutron and positive one at the center. A similar result is obtained with the 
different one-gluon exchange force, acting repulsively between two negatively charged d-quarks. 
But enough about models! 

The starting point is for a model independent analysis is the use of transversely localized 
nucleon states [7J, ||] : 



P(r,n) 



R 2 {r) 



(s\o - r(\ + o - n)o -r\s) . 



(1.3) 
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where |/? + , p, A) are light-cone helicity eigenstates and JY is a normalization factor. The range of 
integration in Eq. (2T) must be restricted to |p| /? + to maintain the interpretation of a nucleon 



moving with well-defined longitudinal momentum [f7|]. Thus we use the infinite momentum frame. 



Using Eq. (|2JJ) sets the transverse center of momentum of a state of total very large momentum 
p + to zero, so that transverse distance b relative to R. can be defined. Thus one defines a useful 
combination of quark-field operators [{7|] : 

O q (x,b) = / d 4-q\ (~,b) q + ( X -,b) e^~, (2.2) 



4?r V 2 

where the subscript + denotes the use of independent quark field operators. The impact parameter 
dependent PDF is defined [Q] as: 

q(x,b) = (p + ,R = 0,X\6 lj (x,b)\p + ,R = 0,X). (2.3) 

This is the basic density that can be obtained, giving the probablitiy of finding a quark at a trans- 
verse position that carries a fraction x of the the protons plus component of momentum. The use of 



Eq. ( |2.l| ) in Eq. ( |2.3| ) allows one to show [g] that q(x,b) is the two-dimensional Fourier transform 



of the GPDH q : 



q(x,b) = J ±*l**H t & =0,x,t = -q 2 ), (2.4) 

One obtains a relation between q(x,b) and the form factor [§] by integrating q(x,b) over all 
values of x. This sets the value of x~ to 0, so that a density appears in the matrix element. If one 
also multiplies by the quark charge e q (in units of e), sums over quark flavors, and uses the sum 
rule relating the GPD to the form factor, the resulting infinite-momentum-frame IMF parton charge 
density in transverse space is 

P (b) = J £e q jd X q(x,b) = f J^F 1 (<2 2 = qV' ,lb . (2.5) 



<i 



We exploit Eq. ( |2.5| ) by using using recent parameterizations [[H], |TI]] of measured form factors 
to determine p(b). The charge densities of the proton and neutron are shown in Fig. 1 of [^j]. The 
surprising feature is the negative central value of the neutron charge density. This results from the 
negative definite nature of F\ [^]. The neutron charge density has interesting features, as shown in 
Fig. 2 of [§] which displays the quantity bp(b). It is the integral of this quantity that integrates to 
0. The neutron charge density is negative at the center, positive in the middle, and again negative 
at the outer edge. The medium-ranged positive charge density is sandwiched by inner and outer 
regions of negative charge. This interesting behavior needs to be better understood. 

One gains information about individual u and d quark densities by invoking charge symmetry 



[ ]12fl so that u,d densities in the proton are the same asd,u densities in the neutron. The results, 



are that the central up quark density is larger than that of the down quark by about 30%. 



Proton Magnetization Density We recently showed [13], that the two-dimensional Fourier trans 



form of the Pauli form factor F2 plays the role of the maagnetization density. 
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3. Measuring the Non-Spherical Shape of the Nucleon 

While the matrix elements of the spin-density operator reveal highly non-spherical densities, 
experimentally determining the proton's non-spherical nature has remained a challenge. Here we 
explain how matrix elements of the spin-dependent density may be measured using their close 



connection with transverse momentum dependent partem densities [14]. 

The field-theoretic version of the spin-dependent charge density operator is a generalization of 
the operator defined in Ref . [Q] : 

p REL (K,n) = J 4^*-*'* vK0)y°(l + y-n 75 m0^; path)^)|^ 0=0 , (3.1) 

where y is a quark field operator and flavor indices are omitted. The quark field operators are eval- 
uated at equal time and accompanied by a path-ordered exponential link operator Jz? ((0, path) 
needed for color-gauge invariance. This matrix element give the probability of finding a quark of 
three momentum K and spin direction n. It depends on the direction of three vectors K,n and the 
direction of the spin polarization. Conpared with the orginally introduced SDD, there is an extra 
7° in front of the term that depends on n. 

Measuring the spin-dependent densities of requires that the system be probed with identi- 
cal initial and final states. But this condition also enters in measurements of both ordinary and 



transverse-momentum-dependent TMD partem distributions. The latter | |15| ] are: 

(3.2) 



<*> [r] (*,K) = J e iK < (P,SM0)rj?(p,S;n-)yr(S)\P,5l) 



where the specific path n_ is that of their Appendix B. The functions <J>I r l depend on the fractional 
momentum x = K + /P + , the trasnverse momentum K, and the proton spin direction. The operator 



r can be any Dirac operator. In particular, the shape of the proton is revealed [14] through 



It is therefore tempting to try to associate an SDD such as that of Eq. ( ]3.1[ ) with TMDs, but one 



difference is essential. Parton density operators Eq. (3.2) depend on quark-field operators defined 
at a fixed light cone time <^ + = <^ 3 + (§ = while our SDD is an equal-time, ^° = 0, correlation 
function. However, a relation between the two sets of operators is obtained |[T4|] by integrating the 
TMD over all values of x setting E, ~~ to zero, and integrating Eq. ( |3~l| ) over all values of K z so that 
i^ 3 = 0. After integration, ^ = for both functions. Computed models for the transverse spin 
dependent densities are shown in Figs. 1,2 of JT4|]. 

The term h\ T causes distinctive experimental signatures in semi-inclusive leptoproduction 



hadron production experiments see the list of references in [|14p. In each of these cases, the mo- 
mentum of the virtual photon and its vector nature provide the analogue of the vector n needed to 
define the spin-dependent density. The hadronic transverse momentum provides the third, Ky. 
Of special interest is the reaction ep f— > e'nX Here the term h\ T causes a distinctive oscillatory 



dependence on the angle 30 — <ps l , [ 16] where <p is the angle between the momentum of the outgoing 
lepton and the reaction plane in the lepton center of mass frame, and fo l denotes the direction of 
polarization with respect to the reaction plane. 
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It is very exciting that experiments planned at Jefferson Laboratory aim to specifically mea- 
sure h\ T [ |l7| ] and therefore determine whether or not the proton is round. We also note that the 
non-spherical shape of the nucleon has been established in lattice QCD by computing appropriate 
moments of impact parameter dependent gpds. See the talk of Zanotti in this workshop. 

Our summary of SDDs is that these are closely related to TMDs. If h\ T is not 0, the proton is 
not round. Experiment can show that the proton is not round. 
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